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We review the Exact Renormalization Group equations of Wegner and 
Houghton in an approximation which permits both numerical and analytical 
studies of nonperturbative renormalization flows. We obtain critical exponents 
numerically and with the local polynomial approximation (LPA), and discuss 
the advantages and shortcomings of these methods, and compare our results 
| with the literature. In particular, convergence of the LPA is discussed in some 

detail. We finally integrate the flows numerically and find a c-function which 
determines these flows to be gradient in this approximation. 



Ph! 1. Introduction 

The Exact Renormalization Group (ERG) is an yet almost unexplored 

approach to non-perturbative computations in quantum field theory (for a detailed 
review of the method see Ref. 4). Recently, some authors have pursued the idea of 
considering the projection of the exact equations on local actions for constant fields.S 
Though restrictive, this approximation mantains a good deal of non-perturbative 
physics - enough to find satisfactory estimates of critical exponent sErQ- and can be 
improved systematically BB 

To start with, we consider one of the physically equivalent formulation of the 
ERG as originally studied by Wegner and Houghton.0 In this contribution we con- 
sider a scalar theory in d dimensions with a single scalar field <p(q) (we work in the 
momentum representation). The ERG equation, describing how the effective action 
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S changes as the high momentum degrees of freedom are integrated out, is 
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- X^ (g) ^^ + ^ +(1 -2-^y/ (g) ^ +const -' (1) 

where r\ is the anomalous dimension of 0, the prime in the first integral above indicates 
integration only over the infinitesimal shell of momenta e~*Ao < q < Ao, and the prime 
in the derivative indicates that it does not act on the 5-functions in dS/d<f)(q). This 
equation is known as the "sharp cut-off" version of the ERG because the integration 
of modes is reduced to a shell. Here we use the approximation proposed in Ref. 5, 
constraining the effective action to have no other derivative pieces than the canonical 
kinetic term, that is, in coordinate space, 

S = J d d x{\{d,<pf + Vm- (2) 

To exploit the ERG Eq. ([I]) within this approximation, Hasenfratz and Hasenfratz 
set up the local, constant mode projection via the action of the operator 

e x aW) (3) 
obtaining the following differential equation for the effective potential 

V(x, t) = -£ ln(l + V"{x, t)) + d- V(x, t) + (1 - - - T))xV\x, t) + const. , (4) 

where A d /2 = [(47r) d / 2 r(ci/2)]- 1 , the dot is a scale derivative d/dt, x is the con- 
stant mode 0(0), and we again refer the reader to the derivation in Ref. 5. In the 
approximation we are using, Eq.(Q) , we actually leave no room for wavefunction 
renormalization, and this turns out to imply that rj = above. For greater ease of 
calculations, we will actually study the equation for f(x,t) = V'(x,t), trivially found 
from the above: 

/(*, t) = Y (i + P(x,t)) + (1 " t] + (1 + t] ' (5) 

with t] already set to 0. The constant Ad can be absorbed by a rescaling of x thus 
disappearing from the equation above, a fact we will make use of later. This is a 
reflection of universality in Eq. (|5]), whereby the shape of /* will depend on but 
the critical exponents will not. 

The numerical treatment of this equation can be carried out in two fashions. 
On one hand, it is possible to use direct numerical solutions to find the fixed-point 
potential. Once the shape of the scale invariant point is obtained, the eigenvalue prob- 
lem of linearized perturbations around it produces the standard critical exponents. 



We shall employ this method in the next section and so obtain the exact solution of 
the differential equation within the approximation we are considering. On the other 
hand, a more analytic approach consists of expanding f(x) as a polynomial series in 
x 2 . This is called the local polynomial approximation (LPA). The complicated differ- 
ential equation for the fixed point reduces then to an algebraic system of equations. 
Solutions are found with less numerical machinery- 

An issue recently touched upon by Morris^ concerns the convergence of the 
LPA to the real solution. This is a subtle problem which, at the moment, hints at 
the necessity of using exact numerical solutions rather than approximations. Whether 
resummation techniques or expansions around_the minima of the potential are appli- 
cable to solve this problem remains to be seenJlim It is also worth mentioning that the 
same authors has managed to extend the approximation @ beyond constant modes, 
that is to order p 2 , and get results that indeed improved upon the zeroth order ones. 

At a more fundamental level, one may consider the long-standing problem of 
irreversibility of renormalization group flows. So far there is no proof of the c-theorem 
is more than two dimensions. The ERG stands as a candidate to progress further on 
this issue and here we have made the smallest attempt in this direction. We have 
taken the polynomial expansion on the projected Wegner-Houghton equation and 
shown irreversibility on its renormalization group flows. This adds a little bit of non- 
perturbative evidence to the validity of the theorem. 

2. Numerical Solutions 

In this section we study the exact renormalization group equations of Wegner 
and Houghton from a purely numerical approach. Fixed point solutions and critical 
exponents found in this way reflect the "best" results these equations will yield, in 
the sense that no ansatz is introduced which could bring in further approximations 
and errors to the final numbers. When we do make the further ansatz of introducing a 
basis expansion (the LPA) for the effective potential in Sec. 3, these results will serve 
as a benchmark to tell us when and how the basis expansion ceases to be a reasonable 
approximation. 

We begin by investigating, in d = 3, the Wilson fixed point solution f*(x) for 
f(x,t) = dV(x,t)/dx. It satisfies the equation 

1 f*" 5 1 

A well-behaved numerical solution of this equation only exists for a specific value 7* of 
the initial condition /*'(0) = 7* (for the other initial condition we take /*(0) = 0). It 
is simple to investigate numerically that if /*'(0) = 7 > 7* then the solution diverges 
at a finite value of x, while if 7 < 7* the solution is unbounded below. Using these 
bad behaviors as a guide we can zero in on 7* with as much precision as our machine 
and software allow us. We find 

7 * = -0.4615413727.. , (7) 

and at the precision above the numerical integration of the equation holds for < 
x < 0.55. This value of 7* agrees to within a 10~ 3 per cent of the value obtained in 
Ref. 5 for the same problem. 



Armed with a numerical interpolating function for f*(x), we can now study 
the linearization of Eq.(^jj) above to determine the critical exponents. The linearized 
equation has a discrete spectrum of eigenvalues. For d = 3 we expect all but one 
of these to be negative, corresponding to the presence of one relevant operator and 
an infinite tower of irrelevant ones, and this is indeed what we find. The numerical 
determination of these critical exponents is again similar in spirit to that of 7* above: 
for eigenvalues slightly off from the correct ones the numerical solutions will either 
diverge up at a finite x, or be unbounded below. As the correct value is approached 
from both sides the well-behaved profile of the eigenfunction is achieved farther and 
farther in x. The first three critical exponents we find are: 

1 

v = — = 0.689458 (8) 

uj = -uj 2 = 0.5953 (9) 
u 3 = -2.84 . (10) 

The last digit in all figures above is uncertain. The first two of the above are 
also in agreement with Ref. 5. It is also in reasonable agreement with calculations 
performed in other methods (field theory calculations, high temperature expansions 
and Monte Carlo methodsS all yield v = 0.630 ± 0.003 and -1.0 < u 2 < -0.5). 

We have also performed the analogous calculations for d = 2.9. There are a 
number of reasons to do this. First, it is only for d < 3 that we expect to find more 
than one non-Gaussian fixed point. At d — 2.9 we indeed find two of them, which is 
what the theory dictates. Secondly, by choosing a dimension close to 3, we can also 
corroborate the results we obtain through an e-expansion in Sec. 3. Finally, within the 
LPA, detecting multiple fixed points for d < 3 involves also discarding a number of 
spurious solutions in a not entirely systematic fashion. Verifying that such spurious 
solutions do not appear here is important to determine that the limitation is not 
in the equation itself, but in the basis expansion, and the solutions we find purely 
numerically can again be used as a guide in identifying the appropriate solutions in 
the basis expansion. 

We find fixed point solutions in d = 2.9 for the following two initial conditions: 

7 * = -0.54598... (11) 
7 6 * = 0.009928... . (12) 

The subindices 4 and 6 indicate the respective fixed points (0 4 and 6 ) . It is worthwhile 
noting that while 74 m 7*, 7g is very small and positive. This is of course a signature 
of a <ft 6 rather than a </> 4 potential profile. 

For this new fixed point we also expect two rather than one relevant operators, 
and a marginally irrelevant one (and infinitely many irrelevant ones). Our e-expansion 
results of Sec. 3 for e = 0.1 give 

wi = 2 (13) 
uj 2 = 1.02 (14) 
cu 3 = -0.2 . (15) 



We have checked the first two of these numerically and have found good agreement. 



3. Local Polynomial Approximation 



In this section we review some properties of exact solutions to equation (|) as 
well as some results on the polynomial approximation for the solutions, the analytic 
approach consisting of expanding f(x,t) as a polynomial series in x. In particular 
we show that reasonable values for critical exponents and correct coefficients of the 
leading order of the e-expansion can be obtained within this approach without much 
machinery. However, arguments will be provided which show that the approximations 
do not converge to the correct limit. Our discussion here is based on results of Refs. 
5,6,9. 

Let us study first properties of t-independent solutions to Eq. (|) with the factor 
Ad/2 being absorbed by a rescaling of x and the function /. As in the previous section 
we choose the initial conditions to be /(0) = and /'(0) = 7. Solutions can be labelled 
by the parameter 7 and in a vicinity of x = can be represented as a series in odd 
powers of x with finite radius of convergence (see below): 

00 

/ 7 (x) = ]T c 2m+1 ( 7 )x 2m+1 , (16) 

m=0 

The coefficients of the expansion obey the following recurrence relations: 

C2m+3{V = -C2m+l(7) ) ( 17 ) 

Qm Qm 

with 

s m = 2(m + l)(2m + 3), q m = 2(m + 1) — md 

and 

g m {c2m+i{l), ■ ■ -ci(7)) = $I c 2m(7)c2(m-o+i(7) [2(m s t . 

1=0 

Using Eq. (|T7D recursively we obtain 03,05,... as functions of c\\ 

m+2 

c 2m+3 (ci) = b i( m + l ) C V ( 18 ) 

1=1 

The coefficients of this polynomial satisfy certain recursive relations, which follow 
from Eq. (pf) and are easy for numerical resolution. However their explicit form is 
rather cumbersome and here we present only the first two of them: 

&i(m+l) = -—him), (19) 
q m 

b 2 (m + l) = —b 2 (m) 

q — m 

1 m 

- — 5>i(0M™-0[2(™-0 + i]*i- (2°) 



Eq. (|T^) can be easily solved and for b\ (0) = 1 gives 



m mi i m 
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where 



2 k 

d ** = ife3l' fc = 2 ' 3 ' 4 ' 



are the upper critical dimensions. 

For large m the asymptotic formula for the coefficients of the expansion ( IB ) 
can be obtained from the relation ([T 



c — w~^^»'"( 1 + (^))- < 22 > 

where the parameter 0(7) is not determined by the leading order terms. This asymp- 
totics shows that the expansion (|l6l) has the finite radius of convergence equal to 

l/^\a(j)\ and is in agreement with the singular behaviour / 7 (x) ~ (2/(x c (j)(d — 



2)))ln(x c (7) — x)p where x c (7) = l/Jafa) is the position of the singularity in the 
complex x-plane closest to the origin. As it has been already discussed in the previous 
section, numerical study of the equation (|J) shows that for certain values 7 = 7* the 
solution / 7 * (x) does not have singularities on the real positive axis and thus is a phys- 
ical fixed point solution. Obviously 7* = is one of such critical values corresponding 
to the Gaussian fixed point solution. It was claimed in Ref. 5 that for d = 4 there are 
no other fixed point solutions. For d = 3 there is a nontrivial solution with 7* given 
by (0) that describes the Wilson fixed point. 

The location of singularities of nontrivial solutions / 7 (x) for d = 3 was analysed 
in Ref. 9. The logarithmic pole x c (7) is real and positive for 7 < 7* and is complex 
but situated close to the real positive axis for 7 > 7*. For 7 approaching the critical 
value 7*, x c (7) moves to infinity and other singularities become dominant. Thus for 
7 = 7* the closest singularities turn out to be at x* = ±r*e ±ld * with r* = 3.12 and 
9* = 0.2577T. 

The representation (|T6| ) suggests that the fixed-point solution to Eq. (||) can be 



found by making the LPA: 

M 

Im{x) = £ c 2m+1 x 2m+1 . (23) 

m=0 



It is clear that the coefficients C2 m +i satisfy the same recurrence relations ([171) or (PL 



for < m < M — 1 plus the additional condition truncating the expansion (]23|) : 

C2M+M) = + l)cf +2 + . . . 61 (M + 1)] d = 0. (24) 

Solutions ci(M) of this algebraic equation through the relations (0) determine the 
coefficients C2 m +3(ci(M)), m = 0, 1, . . . , M — 1, of the polynomial approximation. 



The important features of the LPA can be summarized as follows: 

(1) cm = is always a solution of (p4|) which gives C2 m +i(ci(M)) = and 
reproduces the Gaussian fixed point solution fo(x) = 0. 

(2) From the factorization of the coefficient b\(M + 1), see Eq. (|T9f) , it follows 
that at the upper critical dimensions d = d^ lt , k = 2, 3, . . . cm = is actually 
a double solution, which indicates a branching of fixed-point solutions below these 
critical dimensions. This is in perfect agreement with the multicritical fixed-point 
solutions known to exist below these dimensions. 



(3) The important question is whether the polynomial approximation (|23|) con- 
verges to a non-trivial fixed point solution / 7 *(x) as M — > oo for some sequence of 
Cm satisfying ([24]). The answer seems to be negative and the reason for this is rather 
simple: whereas the true solution / 7 *(rr) possesses singularities, as discussed above, 
the polynomial approximation /m(x) for any M does not. To be more precise, it is 
easy to see that 



M 



\f M (x) - f r (x)\ < l C 2m+l(Cl(M)) - C 2m+ i(7*)||x| 2m+1 + R M (x) 

m=0 

M 

< |ci(M)- 7 *|^|4 m+1 (7*)||x| 2m+1 + J RM(x), (25) 

m=0 

where the residue Rm{x) for \x\ < r* approaches zero as M — > oo. A numerical 
study carried out in Ref. 9 shows that for d = 3 and M > 15 the quantity 5 = 
\ci(M) — 7*| ~ 0.005 and does not decrease with M. This shows that the polynomial 
approximation does not converge to the fixed point solution f 1 *{x). However the 
smallnes of 5 makes the difference to be quite small for x not large. This enables 
us to get an approximation to the solution and, using it, numerical values for the 
first critical exponents (see below) with decent accuracy, while creating an illusion of 
convergence of the method, as claimed in Ref. 6. 

(4) The lower the dimension, the less trustworthy is this approximation or, 
conversely, the larger is the M needed. Altogether, we find that, for any d, some solu- 
tions represent true fixed points while others are spurious. For lower dimensions, the 
number of true nontrivial fixed points increases, but so does the number of spurious 
solutions. There seems to be no wholly systematic way of discarding these spurious 
solutions. 

(5) To further corraborate the usefulness of the polynomial approximation, an 
e-expansion of Eqs. flT8| ) - ( PTf ) about any critical dimension leads to the known e- 
expansion solution. Let us take d = d c ^ %t — e and look for the solution which is of order 
e below the critical dimension d^ lt : 

Cl (M) = ecS 1) (M) + C(e 2 ). 

Since, as it follows from Eq. ([H]), &i(m) = efc^ (m) + 0(e 2 ) for m > k, the solution 
to pfl) for cS X) (M) is given by 



4 1} (M) 



&j 1} (M+l) 
"4 0) (M+ 1)' 



where b 2 (M) is the nonvanishing part of the coefficient b 2 (M) when e — > 0, i.e. 

b 2 (M) = b { 2 ] (M) + 0(e). Using eqs. (g§ and (0) we calculate b { 2 0) and then c?{M). 
We also get from Eq. or flT|) that c 2m+ i(M) = C(e 2 ) for m > k and 

/ M (x) = eK k H 2k ^{x/\ k ) + 0(e 2 ), (26) 

where H 2k -i is the Hermite polynomial and 

K k = (-l) k+1 c?(M)X k 2 k (2k-l)\\, X k = - 2 . (27) 

- 2 

For example, for d^" 1 * = 3 c^M) = 1/20 and does not change with M for M > 3. 
Note that a simple e-expansion of Eq. fl5|) will lead to a linear equation and thus 
cannot furnish the constant K k . At higher orders in e we expect our results not to 
agree with the standard e-expansion since the present approximation (^]) does not 
allow for wavefunction renormalization. 

Having found particular fixed point solutions for some M, we can now study 
how the renormalization flow approaches these solutions by determining the critical 
exponents. To find them, we study small t-dependent departures from some fixed- 
point profile / 7 *(x): 

f(x,t) = f r (x)+g{x,t), (28) 
where again a polynomial ansatz is chosen for g(x,t): 

M 

g(x,t)= ]T ^ift)^ 1 . 

m=l 

When (|28| ) is substituted in Eq. (j5|) and only linear terms in g are kept, we find: 

M 

8 i = Y.Vi3(<?,d)8 j , (29) 
i=i 

where fly is an M x M matrix which depends on the input values q(7*) and d. 

The critical exponents will be given by the eigenvalues of fl. For the 4 fixed 
point we have calculated the critical exponents numerically up to M = 7 for dimen- 
sions between 2 and 4 in steps of 0.1. For 2.9 < d < 4 our results are plotted in Fig. 
1. 

It is worthwhile noting that as d — > 4, the critical exponents merge with the 
tower of canonical dimensions (2, 0, —2, —4, . . .), which are precisely the critical ex- 
ponents of the trivial Gaussian theory at d — 4 (i.e., the canonical dimensions of the 
((f) 2 , (p 4 , 6 , . . .) couplings in d — 4). This is an indication in our setting of the existence 
of a unique (Gaussian) fixed point at d = 4. We furthermore note that for d = 3 our 
two leading exponents 

v = — = 0.656, w 2 = -0.705 

Ll>i 



match fairly well results gotten by other methods (see the previous section). 

It is also possible to perform an e-expansion on the flow equation around a 
critical dimension d = d c k rit . Then, we substitute Eq. fl2"6|) for f y *(x) in Eq.fl28[) and 
make the following ansatz for g(x,t): 

g(x,t) = exp[(^ 0) + euP + t 2 uof ] )t]{g^x) + e 9l {x) + t 2 g 2 {x)). (30) 
For d = d c k nt — e we find that that the critical exponents are equal to 

= ~ e (' ~ 1 - 2{k ~ V-w ) + 0(e2) (31) 

and at the leading order in e 
£=1,2,3,.... 

The LPA also reproduces in the leading order in e the flow solution from the 
vicinity of the Gaussian fixed point fo(x) = to the fixed point Jm(x) given by Eq. 
fl2"Ep. To simplify the formula we choose an initial condition near the Gaussian fixed 
point when all operators except the one, whose critical exponent is proportional to e, 
are zero. Then we get the flow 

1 „ / x 



2^ 



f(x, t) = en k —H 2k _ 1 (—) + 0(e 

1 + ae euJ ^ kt Xk 

where and n k are given by (|27D and u kk = —(k — 1) (see Eq. (|3lD ). The constant 
a is fixed by the initial condition. 

Critical exponents only characterize the flow very close to a particular fixed 
point. Another option offered by the LPA is to study the flow globally by substitut- 
ing Eq. (|I6D , with coefficients c 2m+ i being functions of t now, directly into Eq. (^). 
Matching powers of x in a Taylor expansion leads to coupled nonlinear flow equations 
for Ci(t) in the form: 

c 2 m+i = w 2m +i(c) , i = l to M , (32) 

where the W2 m +i(c) are certain functions of ci, c 3 , . . . c 2m+ i (see Ref. 6) given by the 
differential flow equation (|). Arguably, a polynomial ansatz does introduce a pertur- 
bative element into the essentially nonperturbative nature of renormalization flows 
between distant fixed points, and our approximation very likely misses some features 
of the true flow. However, we believe that, again, the sensible and rich structure that 
emerges does justify the simplification. 

We have solved the nonlinear flow (|32D numerically with M = 3 in d = 3: 

6c 3 



ci = 2ci + 



1 + Ci 



IScl 20c 5 
c 3 = c 3 - — — + 



1 + Ci) 2 1 + ci 
54c| 90c 3 c 5 

C5 = (T+^~(T+^)*- (33) 



The (ci(t), c 3 (t)) subspace of that flow is shown in Fig. 2. 

We note there the presence of a Gaussian and a Wilson fixed point, and a 
unique trajectory leading from the former to the latter. To determine that this flow 
is gradient and permits a c-function description is the object of the next section. 

4. c-Function 

We now study some features of the geometry of the space of local interactions. 
If the beta functions of a theory can be written as a gradient in the space of coupling 
constants, 

/3*(c) = -^£- (34) 

where is a positives-definite metric, we know that the set of renormalization flows 
becomes irreversible. E3 In such a case, there exists a function C of the couplings which 
is monotonically decreasing along the flows: 

dt dc % y dedd> ~ ' y ' 

making their irreversibility apparent, so that recurrent behaviors such as limit cycles 
are forbidden. In two dimensions it is possible to prove that the fixed points of the flow 
are the critical points of C and that the linearized RG generator in a neighborhood 
of a fixed point is symmetric with real eigenvalues (the critical exponents). 

The renormalization group flows found in the previous section are all well- 
behaved. Therefore it becomes natural to ask whether these flows are gradient, i.e., 
whether there exists a globally defined Riemannian metric and a non-singular 
potential C satisfying Eq . (p4j) . The general solution for an arbitrary number of cou- 
plings M would be extremely difficult. However, we find that it is possible to treat the 
case M = 2, namely, the subspace of mass and quartic couplings. The beta functions 
corresponding to the two couplings ci,cs are given in Eq. ( |53"D (where we restrict to 
c 5 = 0). Because of the positivity of c 3 (c 3 is the coefficient of 4 in V and is required 
to be positive for stability of the path integral) it is appropriate to make the following 
coupling constant reparametrization: 

Cl -> m 2 = ci, c 3 -> A 2 = 6A d c 3 . (36) 

In these new variables the beta functions take the form 

dm 2 _ 2 1 A 2 d A _ (4 - d) 3 A 3 

—j—-m +2( 1 + m 2)> -fa- 2 4(l + m 2 ) 2 ( ' 

and the fixed points become 

Gaussian : ( m G' ^g) = (0, 0); 

4-d J2A(A - d) 
Wilson : k, Xw) = {-—^ 1Q _ rf ) • 



Note that the Wilson fixed point merges with the Gaussian one at d = 4, 
similarly to the situation in Sec. 3. Now, by trial and error and considerable guesswork, 
the following solution to Eq.( j34l) can be found: 

C(m 2 , \) = \{l + ™ 2 ) 4 - |(1 + m 2 ) 3 + ^A 2 (l + m 2 ) 2 - (38 ) 

and 

g%3= (1 + m 2 ) ( 4-d ) " (39) 

C(m 2 ,X) has the expected properties of a c-function: i) it has a maximum at the 
Gaussian fixed point, ii) it has a saddle at the Wilson fixed point, and Hi) there 
is only one flow connecting both points ( we have not normalized the c-funtion to 
one for the Gaussian fixed point as often done in the literature). Naturally, this 
description corresponds to our particular parametrization in terms of m and A, which 
implicitly carries a choice of subtraction point. The variation of C between fixed points 
is reparametrization invariant and its positivity amounts to physical irreversibility of 
the flow. A contour plot of C for d = 3 is given in Fig. 3, which depicts the space of 
theories in the basis given by m and A as a hilly landscape. The Gaussian point is at 
the top of the hill (0, 0), whereas the Wilson point lies on the saddle (—1/7, a/24/ 7). 
For the sake of completeness, let us comment that the first mention of irreversibility of 
the renormalization group flow was spelled ouLin the context of perturbation theory 
by Wallace and Zhia.E-a Later, ZamolodchikovtJ proved a theorem in two dimensions, 
the c-theorem, which relates the irreversibility of the flows to the basic assumption 
of unitarity in the Hilbert space of the theory. Several authorsll3 have subsequently 
come to the conclusion that a similar theorem holds in any dimension in perturbation 
theory. More generally, any expansion where the space of theories is reduced to a 
manifold in a space of couplings will accomodate a c-theorem. Our setting in this 
Letter does not clearly fall into this category, due to the appearance of rational 
functions of the couplings in Eq. fl37D , and the explicit construction of the c-function, 
though to first non-trivial order, might be of relevance. 

A systematic approach to the irreversibility of the renormalization group flow in 
the projected Wegner-Houghton equation should rely upon a computation of Zamo- 
lodchikov's metric (i.e. all two-point correlators between composite operators in the 
theory) . This will require an exact renormalization group equation for the generating 
functional equipped with a source for composite scalar fields. 
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Figure Captions 

Fig. 1. Critical exponents for 2.9 < d < 4 corresponding to the relevant, marginal 
and the first two irrelevant operators in the M = 7 approximation. 

Fig. 2. d = 3 Renormalization group flows projected on mass and quartic coupling 
subspace in the M = 3 approximation. c\ is plotted on the x-axis and C3 on the 
y-axis. 



Fig. 3. c-function contour of Eq. (58). The Gaussian point is at the top of the hill 
(0, 0), whereas the Wilson point lies on the saddle (—1/7, \/2A/7). 
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